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ABSTRACT 

A characterization of the optical turbulence vertical distribution and all the main integrated astroclimatic param- 
eters derived from the and the wind speed profiles above Mt. Graham is presented. The statistic includes 
measurements related to 43 nights done with a Generalized Scidar (GS) used in standard configuration with a 
vertical resolution of ~1 km on the whole 20-22 km and with the new technique (HVR-GS) in the first kilometer. 
The latter achieves a resolution of ~ 20-30 m in this region of the atmosphere. Measurements done in different 
periods of the year permit us to provide a seasonal variation analysis of the Cfj . A discretized distribution of the 
typical Cff profiles useful for the Ground Layer Adaptive Optics (GLAO) simulations is provided and a specific 
analysis for the LBT Laser Guide Star system ARGOS case is done including the calculation of the 'gray zones' 
for J, H and K bands. Mt. Graham confirms to be an excellent site with median values of the seeing without 
dome contribution equal to 0.72", the isoplanatic angle equal to 2.5" and the wavefront coherence time equal to 
4.8 msec. We provide a cumulative distribution of the percentage of turbulence developed below H* where H* 
is included in the (0,1 km) range. We find that 50% of the whole turbulence develops in the first 80 m from the 
ground. The turbulence decreasing rate is very similar to what has been observed above Mauna Kea. 
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1. INTRODUCTION 

The Mt. Graham International Observatory (MGIO) is located on Mt. Graham (32°42'05" N, 109°53'31" W), 
Arizona (US) and hosts three telescopes: the Vatican Advanced Technological Telescope (VATT - D = 1.83 m) 
the Heinrich Hertz Submillimiter Telescope (SMT - D = 10 m) and the Large Binocular Telescope (LBT - two D 
= 8.4 m dishes that, when working in an interferometric configuration, can achieve the resolution of a telescope 
with a 23 m pupil size). The study and characterization of the optical turbulence (OT) distribution in space 
and time is fundamental for the ground-based astronomy in the visible up to the near-infrared range to design 
adaptive optics systems and to optimize their performances. The vertical distribution of the OT (i.e. the 
profiles) is the parameter from which all the integrated astroclimatic parameters derive from. 

In this paper we present a study based on measurements of the C% profiles related to 43 nights and obtained 
with a Generalized Scidar (GS) placed at the focus of the VATT on the Mt. Graham summit, around 250 m 
from the LBT. We invite the readers to refer to the paper Masciadri et al. (2010)P for a more extensively and 
complete analysis. In this context, we will try to summarize the most important results obtained and to focus 
on the applications to the Adaptive Optics applications. 

The scientific motivations for such a long-term site testing campaign are: 
(1) To collect an as rich as possible statistical sample of optical turbulence (OT) vertical distribution (Cjy 
profiles) to be compared with simulations obtained with the atmospherical mesoscale model Meso-Nh with the 
aim to validate the model above Mt. Graham. This is a key milestone for the ForOT project whose final goal 
is to predict the optical turbulence above astronomical sites P The measurements from a vertical profiler such 
as a Generalized Scidar are crucial for the validation of such a kind of models. It is our interest to collect a 
heterogeneous sample of measurements taken in different periods of the year and different turbulence conditions 
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Table 1. Classification of the GS campaign measurements. 'GS': C% profiles retrieved from the standard GS. 'WB': C% 
profiles retrieved from the AC frames of the GS measurements obtained with wide-binaries. 'HVR-GS': C% profiles 
retrieved from the wide-binaries auto-correlation (AC) and cross-correlation (CC) frames following the technique described 
in Egner & Masciadri (2007) P 



Sample 


Nights 


Measurements 


Hours 


Resolution 


'GS' 


43 


16657 


163 


AH(0) - 1km 


'WB' 


15 


3659 


6.2 


AH(0) ~ 200 m 


'HVR-GS' 


15 


2812 


5.1 


AH(0) - 25 m 



in a way to better control the model behaviors under different conditions and to better validate the model itself. 
The atmospherical models has been used for the first time to reconstruct and characterize the Cjy profiles by 
Masciadri et alP^ Since that time many progresses have been achieved by our group: the model has been 
applied to different astronomical sites in a simple monomodel configuration^^ aud, more recently, in a grid- 
nesting configuration!^ A new calibration technique has been proposecP and statistically validated^ and the 
first application of the Meso-Nh model as a tool of turbulence characterization has been presented.^ However, 
so far we always could access to GS measurements concentrated in a well defined period of time. The site testing 
campaigns done with a GS at Mt. Graham aimed to comply with the necessity to diversify the measurements 
sample. 

(2) To provide a characterization of all the most important integrated astroclimatic parameters above the site 
of the Large Binocular Telescope (LBT) . 

(3) To provide an as rich as possible statistical sample of the high resolution vertical distribution (AH = 200-250 
m and 25-30 m) of the optical turbulence in the first hundreds of meters up to 1 km to support the feasibility 
studies of new generation instruments for the LBT, such as the LBT Laser Guide Stars system ARGOSp^^ 
that, in its first baseline, is planned to work with a GLACQ configuration. The GLAO efficiency, indeed, strongly 
depends on the turbulence distribution and strength near the ground. The study we intend to perform can be 
achieved thanks to a new technique that has been recently proposed called High- Vertical Resolution Generalized 
Scidar (HVR-G£p2) that aims to reconstruct the optical turbulence vertical distribution in the first kilometer 
above the ground with a resolution 4 up to 10 times higher than what has been done so far with standard vertical 
profilers such as the GS (typically AH ~ 1km). In Egner & Masciadri (2007^ the validity of this technique has 
been proved and it is now our intention to characterize the turbulence distribution in statistical terms, to verify 
how the turbulence decreases in the first hundreds of meters above the ground and provide inputs to test the 
ARGOS performances. 

2. INSTRUMENTS AND DATA-SETS 

Two instruments have been used for this study: the GS and the HVR-GS. We used the GS as developed by 
McKenna et alP^l The HVR-GS technique has been introduced recently^ anc j aims to measure the profiles 
with a high vertical resolution (AH ~ 25-30 m) in the first kilometer above the ground. We refer to Egner & 
Masciadri (2007/^ for the details of the techniques. 

We collected and analyzed observations related to 43 nights (Table 1 in Masciadri et al. (2010) 1 ). Table 
[1] shows the code used to identify the typology of the data-set: 'GS' indicates the standard GS measurements 
extended on ~ 20 km (43 nights), 'WB' indicates the retrieved from the autocorrelation (AC) frames 
associated to wide binaries (AH ~ 200-250 m) (15 nights) and 'HVR-GS' the Cf^ obtained with high vertical 
resolution (AH ~ 25-30 m) (15 nights) in which both the autocorrelation and cross-correlation frames have been 
treated. The samples of the three categories are differently rich because we started to use the HVR-GS more 
recently and this new technique has a higher rate of measurements rejection. We note that, to avoid biases in 
the estimates in the HVR-GS, it has been decided to discard from the statistic all doubtful cases characterized 
by the presence of clouds or cirrus. A method that we called 'normalization^ has been applied to the sample of 
high resolution measurements (15 nights). Among other advantages it permits us to provide a turbulence budget 
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Table 2. Median, first and third quartiles values of the main integrated astroclimatic parameters above Mt. Graham (43 
nights): seeing in the total troposphere, isoplanatic angle, wavefront coherence time, integrated equivalent wind speed. 
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50 th 
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25*' 1 


50 th 


75 th 


25 th 


50 th 


75 th 


e (arcsec) 


0.65 


0.95 


1.34 


0.53 


0.61 


0.72 


0.89 


1.19 


1.50 


9o (arcsec) 


1.6 


2.5 


3.6 


3.1 


3.8 


4.5 


1.4 


2.0 


2.7 


to (msec) 


2.7 


4.8 


8.7 


6.4 


10.1 


14.6 


2.5 


3.8 


6.2 


V (ms- 1 ) 


5.1 


7.2 


9.3 


3.7 


5.1 


7.8 


5.9 


7.7 


9.6 



representative of the whole sample of 43 nights. Basically, the morphology of the turbulence energy distribution 
(shape of the Cjf versus the height) is retrieved from the observation of wide-binaries for 15 nights and with 
the 'normalization' procedure the turbulence energy of the first kilometer detected with the standard GS is 
redistributed in thin turbulent layers according to the profile morphology reconstructed with the CC frames. 

3. INTEGRATED ASTROCLIMATIC PARAMETERS 

Table [5] summarizes the median (50 th ), first (25 th ) and third (75 th ) quartiles for the three main integrated 
astroclimatic parameters calculated for the following groups: the whole sample of 43 nights, the summer and 
the winter time. A composite wind speed profile has been used to calculate the median wavefront coherence 
time To. Below 2 km the wind speed retrieved from the GS has been used; above 2 km the wind speed profile 
retrieved from the European Centre for Medium-Range Weather Forecast (ECMWF) analyses extracted in the 
nearest grid point (32.75°N, 110.00°W) to the Mt. Graham (~ 11.5 km northwest of the summit) has been used. 
The composite method for the wind speed revealed to be the best solution for the calculation of ToP^ We refer 
the readers to Masciadri et al. (2010^ for the calculation of the cumulative distribution of all the integrated 
astroclimatic parameters. Finally, to quantify the contribution of the seeing provided only by the atmosphere, 
the dome seeing (sd), calculated with the method described by Avila et al. (2001)P^has been subtracted from 
the total seeing (e). Knowing that the median seeing in the whole atmosphere (included the dome seeing) is 
e=0.95 arcsec and that ed= 0.52 arcsec, it follows that the median seeing related to the whole atmosphere without 
the dome contribution for the richest statistic we collected so far (43 nights) is £tot= 0.72 arcsec. 

4. COMPOSITE PROFILES FOR ADAPTIVE OPTICS 

The 'composite profiles ' are used to represent efficiently and in statistical terms the vertical distribution of the 
optical turbulence in a discretized number of layers particularly suitable for AO simulations. The method is 
commonly used by many authors in the fielo^ 16 ^ 18 " 20 ! and it consists in identifying a finite number of vertical 
slabs covering the whole troposphere (~ 20 km) with their correspondent value of seeing (Aei or AJ^) so that 
the total turbulence integrated on the troposphere is conserved. J is defined as: 

oo 

J = Jc 2 N (h)dh (1) 



and it is related to the seeing as: 

J = 9-10- 11 -A 1 / 3 -£ 5 / 3 (2) 

where A is expressed in meters, e in arcsec and J in m 1 / 3 . The turbulence in the free atmosphere (h > 1 km) 
and in the boundary layer (h < 1 km) is treated in an independent way to permit to study different combinations 
of probabilities for the OT vertical distribution. Table [3] reports the AJ^ values in the range h > 1 km calculated 
at different heights and obtained from the Cjy associated to the rg related to the 20-30 % of its cumulative 



Table 3. Composite profiles for li > 1 km. In the first column the boundaries of the vertical slabs. In each column is 
reported the value of J that is proportional to the integral of the optical turbulence in the correspondent vertical slab. 
These composite profiles are statistically representative for 43 nights. 
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'Good' 


'Typical' 


'Bad' 
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J (m 1 


/3) 


J (m 1 


/3) 


J (m 1 


/3) 


14000-20000 


7.61e- 


15 


1.02e- 


14 


1.91c- 


•14 


12000-14000 


5.10e- 


•15 


9.65e- 


•15 


1.57e- 


•14 


10000-12000 


7.24e- 


•15 


1.27e- 


14 


2.89e- 


•14 


8000-10000 


8.52e- 


•15 


1.72e- 


14 


3.22e- 


•14 


6000-8000 


6.45e- 


•15 


1.05e- 


•14 


1.77e- 


14 


4000-6000 


9.50e- 


•15 


1.23e- 


•14 


2.10e- 


•14 


3000-4000 


9.18e- 


•15 


1.14e- 


•14 


1.45e- 


•14 


2000-3000 


1.97e- 


14 


3.39e- 


•14 


4.02e- 


•14 


1500-2000 


6.98e- 


•15 


1.52& 


•14 


2.86e- 


•14 


1000-1500 


5.28e- 


15 


1.50e- 


14 


2.86e- 


•14 



Table 4. Composite profiles for h < 1 km after the 'normalization' for the f 3 ; factor. These composite profiles are 
statistically representative for 43 nights. In the last line are reported the J values obtained without the dome contribution 
(median values: £d,25 = 0.35 arcsec, Ed, so = 0.52 arcsec, 6^,75 = 0.70 arcsec). 



Bins 


'Good' 


'Typical' 


'Bad' 


(m) 


J (m 1 / 3 ) 


J ( m V3) 


J (m 1 / 3 ) 


900-1000 


4.35e-15 


7.06e-15 


1.33e-14 


800-900 


2.48e-15 


3.53e-15 


6.29e-15 


700-800 


3.30e-15 


6.04e-15 


1.04e-14 


600-700 


5.72e-15 


1.09e-14 


2.05e-14 


500-600 


4.30e-15 


7.70e-15 


1.38e-14 


400-500 


4.21c-15 


9.60e-15 


1.86e-14 


300-400 


2.05e-14 


4.26e-14 


8.26e-14 


200-300 


6.15e-15 


1.14e-14 


2.53e-14 


100-200 


2.03e-14 


4.14e-14 


8.11c-14 


0-100 


1.79e-13 


3.27e-13 


6.17e-13 


0-100 


3.4e-14 


8.5e-14 


2.20c-13 



distribution ('good' case), to the 45-55 % ('typical' case) and to the 70-80 % ('bad' case). Measurements from 
the sample 'GS' (see Table [TJ are used. 

Again we refer the reader to Masciadri et al. (2010) 1 for precise information on data-reduction analysis. The 
C^r profiles retrieved from the 'WB' sample characterized by a Ah ~ 200 m near the ground (H < 1 km) have 
a suitable vertical resolution to calculate the composite profiles in this vertical range for applications to ARGOS 
(field of view — 4 arcmin) because the turbulence developed below H mi „ = AX/ (2-6) ~ 200 m (AX = 0.5 m is 
the pitch size i.e. the projection of the actuator of the deformable mirror on the pupil of the telescope) is resolved 
by the instrument. Basically we do not need a higher vertical resolution for this application. Table 2] reports 
the composite distribution for h < 1 km equivalent to 43 nights. It includes the dome contribution frequently 
preferable for AO simulations. The median value of the composite profiles (Tableland Table @] central columns) 
permits us to calculate the percentage of the turbulence developed above different heights h with respect to 
the turbulence developed in the whole troposphere (Table [SJleft side) as well as the percentage of turbulence 
developed in the (0, h) range with respect to the turbulence developed in the whole troposphere (Tables-right 
side). Table |6] shows the same calculation obtained in case the median dome seeing is subtracted. 

The composite distribution of the permits also to quantify the extent of the gray zone. We note, firstly, 
that H m ax — t/0 (where r is the residual wavefront coherence size after GLAO correction i.e. the residual 



Table 5. Dome seeing included - Left: Percentage of turbulence developed above the height h with respect to the 
turbulence developed on the whole troposphere. Right: Percentage of turbulence developed between the ground and 
the height h with respect to the turbulence developed in the whole turbulence. Second and fourth columns are obviously 
complementary. 



u J[h,top]/Jtot 

M (%) 

1000 25 

800 26 

600 29 

400 31 

300 38 

200 40 

100 47 



h J[0,h]/J*<rf 

M (%) 

1000 75 

800 74 

600 71 

400 69 

300 62 

200 60 

100 53 



Table 6. Dome seeing subtracted - Left: Percentage of turbulence developed above the height h with respect to the 
turbulence developed on the whole troposphere. Right: Percentage of turbulence developed between the ground and 
the height h with respect to the turbulence developed in the whole turbulence. Second and fourth columns are obviously 
complementary. 
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fwhm_K = 0.13", fwhm_H = 0.2", fwhm_J = 




Figure 1. Extent of the 'gray zone' i.e. H m i„ < h < H mal for different fields of view and wavelengths in the case of median 
distribution (50%). Umax is calculated for J (dotted line), H (dashed line) and K (dot-dashed-dot line) band. H m i„ (bold 
style line) is the same for all the wavelengths. The pupil size is D = 8 m and the pitch size AX = 0.5 m. The FWHM is 
equivalent to r i.e. the residual wavefront coherence size after correction. 



Table 7. Values of H maa; calculated for 9 — 4 arcmin and different residual FWHM. The FWHM values are obtained with 
the GLAO simulations using as inputs Tableland Table 2] H m i„ = AX/20 ~ 200 m for all the wavelengths. 



75% 


FWHM 


^max 


50% 


FWHM 


^max 


25% 


FWHM 






(arcsec) 


(m) 




(arcsec) 


(m) 




(arcsec) 


(m) 
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0.43 


378 
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0.30 


567 
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0.18 


945 
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0.37 


606 
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0.20 


1123 
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0.11 


2042 
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0.25 


1208 
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0.13 


2324 


K 


0.08 


3777 



FWHM) depends on the wavelength. Using the Cff composite profiles extended on 20 km as an input (Table [3] 
and Table 0]), with the combination 'good'- 'good', 'typical'- 'typical' and 'bad'- 'bad', in the boundary layer and in 
the free atmosphere, the ARGOS simulations provide the residual values r (L. Busoni, private communication). 
From these values of r that we reported in Table [7] we can retrieve the values of H m ax for the wavelengths in the 
near- infrared range: J, H and K (Table [7]). Knowing that H m i„ = 200 m, as explained previously, we deduce 
that the 'gray zone ' extends in the (200 m - 378 m) range and assume its smallest value when we observe in J 
band and we consider the 'bad'-'bad' case (75% case). It extends in the (200 m - 3777 m) range and assumes its 
largest value when we observe in K band and we consider the 'good'-'good' case (25% case). Fig. [T] shows the 
H max (blue, red and yellow lines) for different field of view in the case the residual FWHM of GLAO simulations 
has been obtained with the central column of Table [3] and Table [4] i.e (case 50%). 

5. OPTICAL TURBULENCE VERTICAL DISTRIBUTION: C 2 N 

5.1 GS: vertical distribution on the whole troposphere 

Figure [5] shows the median Cff profile obtained with the whole data-set of 43 nights, the summer and the 
winter periods. The morphology of the vertical distribution of the optical turbulence (Cjj profile) shows that 
the greatest turbulence contribution develops in the first kilometer above the ground. Between 1 and 10 km we 
observe a set of minor peaks changing their position and strength during the year. At around 10 km we observe 
the typical secondary Cf^ peak developed at the jet-stream level. For what concerns the seasonal variation we 
observe that the ground layer bump, responsible for most of the turbulence budget, shows a clear seasonal trend 
indicating larger turbulence strength in winter than in the summer period. In the free atmosphere we observe 
the interesting effect of the secondary C% peak located at 10 km in winter time that shifts to higher heights 
(~ 14 km) and is characterized by a weaker strength in summer time. The latter effect (called 'a effect') has 
been put in evidence the first time by Masciadri & Egner (2006)P^ and widely discussed later in Masciadri et al. 
(2010) . We refer the reader to the latter paper for a detailed discussion of the seasonal variation trend in the 
free atmosphere and the 'a effect'. 

5.2 HVR-GS: vertical distribution for h < 1 km 

The optical turbulence vertical distribution with high resolution (20-30 m) in the first kilometer is obtained 
with the method called HVR-GS presented in Egner & Masciadri (2007) PI The HVR-GS data-set is obtained 
taking the integral of the profiles retrieved from the AC frames and redistributing the energy in the first 
kilometer according to the detected triplets in the CC frames. Three different strategies can be used to study the 
turbulence spatial distribution in the boundary layer. The usefulness of each method depends on the application 
one intends to give to the analysis. We study: 

(A) the median of the (and/or J) profiles. 

(B) the average of the Cfj (and/or J) profiles. This method is very useful for comparisons of measurements and 
simulations obtained with atmospherical models. Moreover this operator has the advantage that the mean of Ji 
is equal to the J retrieved from the mean of the C^ i (where Ji and C^ i refer to the each individual profile). It 
is not the case for the median. 

(C) the composite profiles as calculated in Section @] This method is very useful for applications to Adaptive 
Optics. 





Figure 2. Median C N profile obtained with the complete sample of 43 nights, the summer [April- June] and winter 
[October-March] time samples. Results are obtained with the standard GS technique. 



Following the strategy (A) we find that C% = for h > 25-30 m. This can be explained with the fact that 
each vector (C^- profile) has many zeros. This is due to the fact that, during the monitoring, the system detects 
spikes and thin layers that changes position and duration with the time. In other words, the use of the median in 
these cases, can be misleading and it can provide very low values of the C]^ . Considering that the strategy (A) 
is not really useful to characterize these measurements in our case, the strategies (B) and (C) have been used. 
Figure [3J- left shows the result obtained following the strategy (B) i.e. the mean of the Cfj profiles (Ah ~ 25-30 
m) calculated from the sample 'HVR-GS'. It is therefore representative of 43 nights. The two profiles (with and 
without dome contribution) are shown in proximity of the ground (Figure [3J- right). To retrieve the typical scale 
height B of the exponential decay of the mean Cjy profile the measurements done below 125 m have been fitted 
with an exponential law (Eq.Q) as we already did in Stoesz et al. (2008) 



y = A ■ e<-VB) (3) 

where A and B are free parameters. The calculation is obviously shown only in the case in which the dome 
contribution is subtracted. The fit gives A = 3.34-10 -15 and B = 37.4 m (Fig J3J- left). If we limit the analytical 
fit to the first 30 meters, the scale height B = 28 m. The exponential decay is however an absolutely arbitrary 
analytical law. The important issue is that these results definitely indicate that the HVR-GS technique is able 
to put in evidence that the turbulence decays above typical astronomical sites in stable night time conditions in 
a much sharper way than what has been predicted and quantified in the past. Indeed, the Hufnagel model in 
proximity of the surface^ in night conditions states that the turbulence scales as h~ 2 / 3 (FigJH]). For the Hufnagel 
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Figure 3. Left: Mean C% profile (A/i~20-30 m) calculated with the sample 'HVR-GS' and the dome seeing contribution 
included (black line). Hufnagel model (red line). Right: Zoom of the first 200 m. Near the ground it is put in evidence 
the different C% shape in case the dome contribution is included (black line) and excluded (green line). 



model (Fig|3]-left) the Cjy decreases of one order of magnitude within 1 km while our results indicate that the 
C^f decreases of one order of magnitude within ~ 60-70 m (FigJSJ-right). It remains interesting the (800 m - 1 
km) range in which visibly a weak turbulence develops. One should expect a smoother connection between the 
boundary layer and the free atmosphere. We think that this is just an artifact effect due to the use of a quite 
different resolution below and above 1 km. From the point of view of the AO simulations this small gap should 
not cause any problems. It should be enough to implement a weak convolution to slightly smooth out the Cff 
vertical profiles at the interface located at 1 km so to obtain a less abrupt connection of the profile above 
and below 1 km. 

We note that the presence of no signal (therefore C% = 0) might potentially indicate that the turbulence 
strength is not equal to zero but simply weaker than the threshold Cf^ =10~ 16 (associated to an equivalent 
J=2.5T0 -15 ). In Masciadri et al. (2009)p2] ft uas been calculated, in a post-processing phase, the most con- 
servative case in which we assigned C^- =10~ 16 where there is no signal. After a more careful investigation we 
observed that the latest distribution is associated to a too large total J in the boundary layer (the turbulent 
energy should not be conserved) and, for this reason, it can be discarded. 

Figure H] shows the results obtained following the strategy (C) that is the mean C]^ profiles associated to the 
J (or r ) related to the 20-30 %, 45-55 % and 70-80 % of the cumulative distribution. In Masciadri et al. (2010^ 
the numerical values for the correspondent J values are reported. Curiously the first grid point near the ground 
of the 75% case distribution shows a weaker value with respect to the 25% and 50% cases. This is due to the 
fact that the third quartile dome seeing (£^75), that has been subtracted from the original value is particularly 
large (0.70 arcsec). The morphology of the C% in the first kilometer is very interesting showing several thin 
layers and a very weak turbulence between 800 m and 1 km. We note that the sample on which we calculate the 
average in each slot [20-30]%, [45-55]% and [70-80]% is of a few hundreds of Cff profiles therefore there are no 
doubts that this structure reproduces some real distribution. In terms of morphology of the turbulence profile we 
find therefore that the higher the vertical resolution the thinner the size of the detected layers. This conclusion 
is perfectly coherent with the turbulence structure resolved by balloons equipped for the Cjf measurements.^ 

Figure [5] shows that the percentage of turbulence P(h) developed in the (0, h) range with h in the (0, 1 km) 
vertical slab. The function P(h) is defined as: 
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Figure 4. Mean C% profiles calculated from the corresponding J (or ro) values related to the 20-30 %, 45-55 % and 70-80 
% ranges of the J cumulative distribution. The dome contribution is excluded. In Annex C the numerical values for the 
correspondent J values are reported. 
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(%) (%) 

Figure 5. Percentage of turbulence developed between the ground and the height h with respect to the turbulence developed 
in the whole atmosphere (~ 20 km) as retrieved from the HVR-GS measurements and extended to the first kilometer. 
On the left hand is shown a zoom of the picture centred on the first hundreds of meters. 



where the C% profile is that associated to the 45-55 % case (Fig. HJ-centre). It is worth to note that the error 
bars for the HVR-GS technique is of the same order of half of the vertical resolution i.e. ± 12-15 m. This derives 
mainly by the definition of the zero point i.e. the ground that is characterized by an uncertainty equivalent to 
half of the vertical resolution. 

From FigJSJwe retrieve that around 50% of the turbulence developed in the whole 20 km is concentrated below 
80±15m. This result is substantially different from the preliminary indications obtained in Egner & Masciadri 
(2007^ and the turbulence seems to be much more concentrated near the ground. We also note that the Cf^ 
morphology below 1 km decreases in a different way if we look at the sample 'WB' and the 'HVR-GS'. This is 
not a contradiction and it can be explained with the fact that the turbulence spatial distribution for the 'WB' 
and the 'HVR-GS' samples is not necessarily the same because the first is based on the AC frames while the 



second is based on the CC frames. It is interesting to note that the higher the vertical resolution, the sharper is 
the exponential decay of the morphologic turbulence structure. 

We note that, simultaneously to our study, some other authors^ recently investigated the turbulence structure 
near the ground at high vertical resolution. Even if the instruments employed were different these studies present 
many similarities in the results. From a qualitative point of view, we note that, also in that case, the turbulence 
appears well confined near the surface. From a quantitative point of view, things are more delicate. Looking at 
Table 3 (in that paper) it appears that their data-reduction is more similar to our method (C) than to the others 
methods. Figure H2 shows the J-profile retrieved from Fig. 0](5O% case) overlapped to the J-profile retrieved from 
Table 3-Chun et al. (2009). ^ We note that above Mt. Graham the turbulence vertical sampling is 25 m, above 
Mauna Kea the sampling increases from 15 m up to 80 m and it extends only up to 650 m. The turbulence 
vertical distribution appears very similar. Above Mauna Kea a local minimum is present at ~ 45 m above the 
ground more or less in correspondence of the abrupt detection break due to the sensitivity threshold from LOLAS 
(Table 3-Chun et al. (2009)% We highlight that the evident huge 'turbulent vacuum zone ' between 560 m and 
1 km (SLODAR) simply means that the turbulence is not measured in this vertical slab, not that turbulence is 
not present. 

6. CONCLUSIONS 

In this paper we present the results of a study aiming to characterize the optical turbulence at Mt. Graham. We 
present a general overview of the statistics (43 nights) of the Cjj profiles and all the main integrated astroclimatic 
parameters and their seasonal trends. The main conclusions we achieved are: 

(1) With a median seeing e= 0.95 arcsec (e= 0.72 arcsec without dome contribution), isoplanatic angle 6*o= 2.5 
arcsec and a wavefront coherence time to=4.8 msec, Mt. Graham confirms its good quality in terms of turbulence 
characteristics typical of the best astronomical sites in the world. All the integrated astroclimatic parameters 
(the seeing, the isoplanatic angle, the wavefront coherence time and the equivalent wind speed) show a clear 
seasonal trend that indicates better turbulence conditions and weaker equivalent wind speed Vo in summer with 
respect to the winter. 

(2) The ground layer is characterized for the first time with a high resolution (200-250 m and 20-30 m). The 
turbulence exponentially decays above Mt. Graham with a much sharper profile than what has been supposed 
so far and expressed with the Hufnagel model. Three different strategies of analysis aiming to investigate 
the morphology of the turbulence spatial distribution have been presented. We find that around 50% of the 
turbulence developed in the whole atmosphere is concentrated below 80±15 m from the ground and 60% of the 
turbulence in the first kilometer. This evidence together with the favorable large 9$ observed above Mt. Graham 
(particularly in the spring/summer time) represent extremely favorable conditions for astronomical observations 
assisted by a LGS/GLAO system such as ARGOS. 

(3) We observe that the higher is the vertical resolution of the tool used to measured the turbulence vertical 
distribution the sharper is the turbulence decreasing. 

(4) The percentage of turbulence developed below the primary mirror of the LBT (i.e. ~ 35 m from the ground) 
is around 33%. However this estimate has to be considered with precaution because the uncertainty (2cr~ 25-30 
m) is of the same order of magnitude of the vertical resolution (AH ~ 25-30 m) and in the first hundred meters 
the turbulence decreases very sharply. 

(5) It appears evident that at Mt. Graham the turbulence decreases above the ground similarly to what 
observed above Mauna Kea in a more or less simultaneous study (Chun et al. 2009) 2 "-" performed with different 
instrumentation. 

(6) A composite distribution of the turbulence on the whole 20 km is calculated to be used as input of AO 
simulations of the LBT Laser Guide Star system named ARGOS and the calculation of the gray zones for the 
near-infrared J, H and K band is done. The gray zone extends from a minimum of (200m - 378 m) in J band 
with the 'bad '-'bad' case up to a maximum of (200m - 3777 m) in K band in the case 'good'-'good' case. 

(7) A clear seasonal variation trend has been observed in proximity of the ground and in the jet-stream 
regions. These measurements confirm the first evidence of the seasonal trend observed by Masciadri & 
Egner (2006)^ above other astronomical site. The physical model proposed by Masciadri & Egner (2006)^ that 
is able to explain the seasonal effect of the secondary peak of the Cfj called 'a effect', is confirmed and refined. 
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Figure 6. Red line: turbulence vertical distribution (J-profile, 50% case) above Mt. Graham as calculated in Fig|4]-centre. 
Black line: turbulence vertical distribution (J-profile, 50% case) above Mauna Kea as calculated by Chun et al. (2009) in 
Table 6-centre. Mauna Kea measurements extend up to 650 m. 



(8) For the first time we observed a seasonal trend of the dome seeing. This is certainly a topic that deserves a 
more careful investigation in the future. 
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